Unzipping flux lines from extended defects in type-II superconductors. 
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With magnetic force microscopy in mind, we study the unbinding transition of individual flux 
lines from extended defects like columnar pins and twin planes in type II superconductors. In the 
presence of point disorder, the transition is universal with an exponent which depends only on 
the dimensionality of the extended defect. We also consider the unbinding transition of a single 
vortex line from a twin plane occupied by other vortices. We show that the critical properties of 
this transition depend strongly on the Luttinger liquid parameter which describes the long distance 
physics of the two-dimensional flux line array. 



The competition between thermal fluctuations, pin- 
ning and interactions of vortices in type-II high- 
temperature superconductors leads to many interesting 
physical phenomena . These include the melting of the 
Abrikosov flux-lattice into an entangled vortex-liquid Q 
and the proposed existence of low temperature Bose-glass 
3] , vortex glass .4] and Bragg glass [f| phases. 

Experimental probes range from decoration, transport 
and magnetization measurements, neutron scattering, 
electron microscopy, electron holography to Hall probe 
microscopes. While these experiments yield a wealth 
of information, the possibility of manipulating a sin- 
gle vortex, for example using magnetic force microscopy 
(MFM) 0, has received less attention. Such experiments 
could give rise to a direct measurement of microscopic 
vortex properties which, up to now, have been under 
debate or assumed. The possibility of performing such 
experiments is akin to single molecule experiments on 
motor proteins, DNA, and RNA which have transformed 
biophysics during the past decade by opening a window 
on phenomena inaccessible via traditional bulk experi- 
ments Q. 

Such an experiment was recently proposed by Olson- 
Reichhardt and Hastings These authors suggested 
using MFM to wind two vortices around each other, thus 
probing directly the energetic barrier for two vortices to 
cut through each other. A high barrier for flux line cross- 
ing has important consequences for the dynamics of the 
entangled vortex phase. 

In this Letter we analyze theoretically several possible 
experiments in which a single vortex is depinned from 
extended defects using, for example, MFM. We first con- 
sider the situation where MFM is used to pull an iso- 
lated vortex bound to common extended defects such as 
a columnar pin or a twin plane in the presence of point 
disorder. Using a scaling argument we derive the dis- 
placement of the vortex as a function of the force exerted 
by the MFM near the depinning transition in an arbitrary 
dimension d. We argue that the transition can be char- 
acterized by a universal critical exponent, which depends 
only on the dimensionality of the defect. Hence unzipping 
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FIG. 1: The experimental setup considered in the text with 
a columnar defect in 2+1 dimensions. The MFM is applying 
a force in the x direction. 



experiments from a twin plane directly measure the free- 
energy fluctuations of a vortex in the presence of point 
disorder in d = 1 + 1 dimensions. To the best of our 
knowledge, there is only one, indirect, measurement of 
this important quantity in Q. Our conclusions are sup- 
ported by numerical and analytical calculations. Related 
results apply when a tilted magnetic field 10] is used to 
tear away vortex lines in the presence of point disorder. 

Next, we study the effect of interactions between vor- 
tices on the unzipping of a single vortex. We consider 
a system in which vortices are preferentially bound to a 
thin two dimensional slab, such as a twin plane, of say 
YBCO 1], in a three dimensional sample. A similar sit- 
uation can be achieved by artificially inserting (with, for 
example, molecular beam epitaxy) a thin plane with a 
reduced lower critical magnetic field H c i into a bulk su- 
perconductor; we require that the density of vortices in 
the twin plane is much higher than in the bulk of the 
sample. We imagine that MFM is used to pull a single 
vortex out of the slab, thus creating an effective mag- 
netic monopole inside it. The displacement of the vor- 
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tex from the slab as the transverse force exerted by the 
MFM is increased then depends on the physics of the 
two dimensional vortex liquid which resides in the slab. 
Specifically, the "Luttinger liquid parameter" which 
controls large-distance behavior of the vortex liquid also 
controls critical properties of the unbinding transition of 
the vortex from the slab. We study the "unzipping of 
Luttinger liquids" both with and without point disorder, 
and argue that this setup can be used as a sensitive probe 
of the two dimensional physics. 

Consider first the unzipping of a single vortex from an 
extended defect. We describe the flux line by a function 
r(r), where r denotes a coordinate parallel to the defect 
(see Fig. H), and r denotes the (d — l)-dimensional co- 
ordinates perpendicular to it. In the absence of external 
force the appropriate free energy for a given configuration 
r(r) of the vortex is given by |l|: 
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Here 7 is the line tension and L is the length of the 
sample along the r direction. The point disorder /x(r, r) 
is assumed uncorrelated and Gaussian distributed with 
/I = and /i(r, r)/i(r', r') = cr6(v — t')S(t— r'), where the 
overbar denotes an average over realizations of the dis- 
order. V(r) is a potential describing the d'-dimensional 
extended defect. 

To study unzipping, we consider a total free energy 
T = J-q + Fi, where T\ = — f • r(L) = — / Q L f • d T v(r) dr. 
Here f stands for the local force exerted by the MFM in 
the transverse direction and we assume that r(0) = 0: 
below the unzipping transition the flux line is unaffected 
by the boundary conditions at the far end of the sample. 

As the force, f , increases the free energy density of the 
unzipped portion of the vortex decreases. In contrast, 
the free energy density of the bound part is, clearly, in- 
dependent of the force. The vortex will be unzipped when 
/ = f c such that the free-energy densities of the bound 
and unzipped states intersect. We study the behavior 
of s = (m), where m is the length along the r direc- 
tion which is unbound as the force approaches the crit- 
ical value f c . The angular brackets and overbar denote 
thermal and disorder averages, respectively. Although 
we focus on disorder averages, particular disorder config- 
urations will be characterized by a pattern of jumps and 
pauses in the unzipping |l2|. Disorder averages can be 
probed experimentally by unzipping many different vor- 
tices. It is straightforward to show that (m) is related 
linearly to the transverse displacement of the vortex at 
the top of the sample, as measured by MFM. We take the 
force to always act in the x direction. Without disorder 
it is known that s ~ (f c — in any dimension and for 
any dimensionality of the defect • 

The universal properties of the unzipping transition 
with point disorder can be analyzed using a simple scaling 



argument adapted from for the unzipping of DNA. 
Consider the free energy J-(m), with contributions from 
three sources. The first is linear in m and is due to the 
average free energy difference between a vortex on the de- 
fect and in the bulk of the sample (/ c — f)m. In addition 
there is a contribution, STd', from optimized free-energy 
fluctuations of the segment of the vortex which is bound 
to the defect. This term arises from point disorder which 
is localized on or near the defect and is expected to be- 
have as m 1 / 2 at large m for a d' = 1 dimensional defect. 
For d' > 1 we expect SJ-41 cx m w ^ d \ where uj(d!) is the 
exponent associated with the free-energy fluctuations of 
a directed path in the presence of point disorder in d' 
dimensions |13| . Finally, there is the interaction of the 
unzipped vortex with bulk point disorder, SJ-^ oc m"^, 
where d > d! is the dimensionality of the sample. Col- 
lecting these terms gives: 



T(m) = a(f c - f)m - bm"^ 



cm 



'(d) 
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Here a, b and c are positive constants and the negative 
signs have been chosen since the behavior is expected to 
be dominated by the minima of the two random poten- 
tials. The exponent uj(d) has been studied extensively in 
the past and it is well known that u(d') < u>{d) for any 
d! < d 13] . For example u>(d' = 2) = 1/3 (twin planes) 
and u(d = 3) ~ 0.22 (bulk sample). Therefore, disorder 
on or close to the defect controls the unbinding transition 
for any dimension and the problem is equivalent to unzip- 
ping from a sample with disorder localized on the defect. 
In practice, disorder is likely to be concentrated within 
real twin planes and near columnar damage tracks cre- 
ated by heavy ion irradiation, strengthening even more 
the conclusions of this simple argument. By minimiz- 
ing the free energy with respect to m we find the typical 
value of length of the unzipped part of the flux line, which 
determines s — (m): 
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v= [l-w(d')]- 1 



(3) 
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Thus we get v — 2 (in agreement with f° r a colum- 
nar pin and a new result, v = 3/2 for a twin plane. 

The heuristic argument leading to Eq. J3J) can be tested 
in a number of ways. First consider unzipping from an 
attractive hard wall in (1 + l)-dimensions with point dis- 
order in the bulk. Upon applying an imaginary-gauge- 
transformation flfl ] (possible below the unzipping tran- 
sition) to the Bethe-ansatz solution for / = 0, obtained 
by Kardar (see Eq. (3.12) in Ref. H3)j we obtain 



dy ny 



exp[-(A- f/j- n)y\ 
[1 - exp(-Ky)] 



(4) 



Here k = a/2 and A is the inverse localization length 
characterizing the pinning to the wall in the absence of 
point disorder. For n <C A — f /j, Eq. reduces to the 
clean result: s ~ (A — //7) _1 , while in the opposite limit 
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one has s ~ (A — k — f /j)~ 2 as expected from Eq. |3] 
This result confirms our expectation that, even though 
disorder is present everywhere, the universal properties of 
the unzipping transition are determined by an effective 
disorder potential generated in the vicinity of the wall 
with v = 2. 

To check the case of a symmetric attractive potential, 
representing a columnar pin, we performed numerical 
simulations for a lattice model in d = 1 + 1. The re- 
stricted partition function of this model, Z(x,t), which 
sums over the weights of all path leading to x, t, satisfies 
the recursion relation 

Z(x,r + l) = 4,o(e y -l)Z(0,r) 

+e n(x,r+i) [ Je f Z (x-l,T) + Je- f Z(x+l,T)] .(5) 

Here fi(x, r) was taken to be distributed uniformly in 
[— Uq,Uq] with a variance a = Uq/3. The variable J 
controls the line tension (J = exp7/2) and was set to 
be J = 0.1 while we used V = 0.1. We work in units 
such that ksT = 1, where ks is the Boltzmann con- 
stant. The partition function was evaluated for each vari- 
ance of the disorder for several systems of finite width 
w = 2L X averaging over the time-like direction (typi- 
cally t ~ 10 6 "time" steps) with the initial condition 
Z(0, 0) = 1 and Z(x, 0) = for x ^ 0. In the vicinity of 
the transition we expect the finite size scaling prediction 
s = L x f [L x (f c — /)"]. For finite L x we expect a smooth 
interpolation between v = 1 (the clean result) and v = 2 
(the asymptotic result with disorder) with increasing L x 
or increasing strength of disorder. Using standard meth- 
ods of finite-size scaling we extracted the exponent for a 
given strength of the disorder as a function of L x . The 
exponent for each value of L x is obtained from the best 
collapse of the data of two systems sizes L x and L x /2. As 
shown in Fig. [5J the data is consistent with v saturating 
at v — 2 for large systems. The crossover to v = 2 is 
much more rapid if the point disorder is enhanced near 
the columnar pin (see the inset in Fig.^J, as might be ex- 
pected for damage tracks created by heavy ion radiation. 

Numerics in d = 2 + 1 dimensions with a twin plane 
or a columnar pin proved more difficult. However, given 
that the universal properties of the unbinding transition 
are determined only by the disorder on the defect, we 
have performed numerical simulations of the model with 
no bulk disorder and ignoring excursions of the vortex 
into the bulk for t values where it is pinned to the de- 
fect. This allows very large systems to be studied. The 
d = 1 + 1 problem, with no bulk disorder, was studied in 
the context of DNA unzipping which yields s oc (f c —f)~ 2 
in agreement with Eq. [3] 12]. For additional analytic re- 
sults, using the replica trick see |lfil |. In Fig. we show a 
data collapse with the anticipated exponent v — 3/2 for 
unzipping from a disordered twin plane [d = 2 + 1), using 
the finite size scaling form s = Lg[L(f c — f) u ], with L 
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FIG. 2: Effective exponent 1/u versus L x for a fixed strength 
of point disorder a = 0.03. The results are consistent with 
the general argument that this exponent should saturate at 
v — 2 as L x — > <x. The inset shows the same exponent vs a c , 
the variance of additional point disorder placed directly on 
the columnar pin extracted from two system sizes L x — 600 
and L x = 1200. It is clear that v — > 2 as a c increases. 
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FIG. 3: Data collapse for unzipping from a disordered twin- 
plane with v = 3/2. See the text for more details. 

the size of the system in the r direction. For the sim- 
ulations we used a model like that in Eq. (JjJJ with an 
in-plane hopping amplitude J = 1/2 and a strength of 
point disorder, a = 4/3. To obtain Fig. EJwe evaluated 
numerically the free energy of the bound part of the vor- 
tex line for each realization of point disorder and used an 
analytical expression for the free energy of the unbound 
part (F u (m)). The latter is readily obtained from Eq. ||TJ 
by summing over weights of paths of a directed random 
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walk which is pulled by a force /: F u (m) = — f 2 m/2^. 
This allows a very efficient calculation of s including av- 
erages over point disorder. The results agree very well 
with the prediction of Eq. ©. 

Consider now unzipping from a two dimensional plane 
with many vortices in a three dimensional sample which 
is essentially free of vortices (see Fig. 0J. A MFM tip 
pulls the top end of one of the vortices out of the plane 
with a force f = fx. As in the single vortex problem, we 
expect an unzipping transition for f > f c such that the 
vortex is pulled out of the two dimensional slab. We again 




FIG. 4: Schematic representation of an MFM tip pulling out 
a vortex out of a 1+1 dimensional flux line array. 

write the free energy of the vortex as a sum of two con- 
tributions: F u (m), which arises from the vortex segment 
separated from the two dimensional slab, and Fo(m), the 
free energy change of the vortices which remain inside. 
Here m is again the length along the r direction which is 
unbound from the plane. 

The free energy F u (m) is evaluated analytically as dis- 
cussed above. As for F (m), clearly, there is a linear con- 
tribution associated with the length to removed from the 
attractive potential of the slab. In addition, a non-linear 
contribution in r arises due to the dislocation created in 
the vortex array (see Fig. 0J. The energy of this dislo- 
cation is determined by the long distance properties of 
a two-dimensional vortex liquid. These are known to be 
controlled by a single parameter g, which is related to the 
elastic moduli of the vortex lattice and temperature |ll| . 
For example, vortex density oscillations without point 
disorder decay as a power-law in distance with an expo- 
nent r] = 2g. Due to an analogy with one dimensional 
quantum bosons, g is often referred to as the "Luttinger 
liquid parameter" . Assuming that vortex lines must exit 
normal to the interface, the energetic cost of the disloca- 
tion, Td-, can be calculated using the method of images: 
Td = ^ff In to, for large to [lr| . 

Having obtained Fo(to) and F u (m) it is straightfor- 
ward to calculate (m) = s, (related to the mean dis- 



placement of the MFM tip by x m = fs/j) as a function 
of / using the free energy 

F(m) = r(f c - f)m + ^1 ln(m). (6) 

where r is a positive constant. At the transition point, 
/ = f c and e~ F ( m ^ T is a power law in to, so the results 
are sensitive to the value of g. For g > 1/4 we find 
the same divergence as for a clean twin plane without 
interactions, 
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(7) 



In contrast, for 1/8 < g < 1/4 one finds a continuously 
variable exponent governing the transition 
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Finally, for g < 1/8 we find that s remains finite as 
f —* fc an d the unzipping is presumably discontinuous. 
Higher moments of to will also be sensitive to the value 

of g- 

Finally, we consider the case of unzipping from a plane 
with many vortices in a three dimensional sample in the 
presence of point disorder. Here we expect the results 
to be sensitive to the boundary conditions imposed on 
the plane. When the boundary conditions are such that 
the number of flux lines is conserved (e.g. in a cylindri- 
cal geometry) we expect the free energy fluctuations to 
behave as 8T{m) oc to 1 / 2 due to the highly constrained, 
almost one dimensional, behavior of the flux lines. Thus, 
the universal exponent is expected to be v — 2 for any 
value of g. In slab geometries however, where flux lines 
can freely enter from the sides of the plane, the single 
flux line contribution to the free-energy fluctuations can 
be screened and the effective fluctuations come only from 
the dislocation energy. Near the vortex glass transition 
at 9 = 1> [4I Il5l lri| the calculation can be carried out us- 
ing the replica approach. Although a rigorous treatment 
requires the method of images with reflected disorder, 
we do not expect this to affect the results and performed 
the calculations with uncorrelated disorder. We find that 
while the quenched average of the free energy is indepen- 
dent of the disorder, its fluctuations grow as (ln(TO)) 1 / 2 
for g > 1 and (1 — g) ln(m) for g < 1. Thus, in contrast 
to the single vortex problem and to the case when the 
number of vortices in the sample is conserved, here fluc- 
tuations in the bulk of the sample are dominant, growing 
as toM 3 ). Using w(3) ~ 0.22, this lead one to the predic- 
tion s ~ (f c — f) -1 ' 28 , allowing a direct measurement of 
w(3). 
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